Reconstructed Rough Phases During Surface Growth 
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Flat surface phases are unstable during growth and known to become rough. This does not ex- 
clude the possibility that surface reconstruction order persists in rough growing surfaces, in analogy 
with so-called equilibrium reconstructed rough phases. We investigate this in the context of KPZ 
type dynamics, using the restricted solid on solid model with negative mono-atomic step energies. 
Long range reconstruction order is strictly speaking absent in the thermodynamic limit, but the 
reconstruction domain walls become trapped at surface ridge lines, and the reconstruction order 
parameter fluctuates critically with the KPZ dynamic exponent at finite but large length scales. 
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Equilibrium surface phase transitions have been stud- 
ied in great detail during recent decades. Various types of 
critical behaviors emerged both in theoretical models and 
experimental systems Dynamic non-equilibrium as- 
pects are still less well understood, in particular whether 
any of those equilibrium transitions persists in the sta- 
tionary state properties of growing interfaces. For ex- 
ample, equilibrium crystal surfaces undergo roughening 
transitions from macroscopic flat to rough structures, 
while growing surfaces are believed to be rough under 
all circumstances |^,^, as confirmed by numerous stud- 
ies of, e.g., KPZ type dynamics Q and other dynamic 
universality classes. Still, it is custom to identify dis- 
tinct growth regimes within this rough phase. So-called 
step-flow type layer-by-layer growth is an example. At 
low temperatures, well below where the equilibrium 
surface roughens, a new terrace is nucleated with an ex- 
ponential small probability. The time scale at which 
this nucleus grows by particle adhesion at its edge into a 
macroscopic domain is much shorter than the time scale 
at which a nucleus for the next layer appears on top of 
this terrace. So at small enough length scales the surface 
looks flat and seems to grow layer-by-layer. In this letter 
we introduce a similar transient phenomenon for surface 
reconstruction order in growing rough phases. 

Surface reconstruction is usually associated with flat 
interfaces. However, roughness not necessarily destroys 
the reconstruction order. Equilibrium reconstructed 
rough (RR) phases are known to exist for so-called mis- 
placement type reconstruction [^J^. The compatibility 
of surface reconstruction with surface roughness depends 
on intricate topological aspects. For example, in missing 
row reconstructed FCC (110) facets the reconstruction 
couples strongly to the surface roughness such that sur- 
face roughening simultaneously destroys the reconstruc- 
tion order For other symmetries, like simple cubic 
(SC) (110) missing row reconstruction, they decouple. 

The RR order parameter must be formulated with care 
1^. In flat SC missing row structures for example, the re- 
construction order can be formulated in two ways, which 
seem equivalent at first, but at the roughening transition 



only one of them vanishes. One formulation keeps track 
of whether the even or odd rows are on top. The other 
one measures it by the (striped) antiferromagnetic or- 
dering of the parity type Ising variables Sr = exp(«7r/ir), 
with hr the surface height. Step excitations in this sur- 
face belong to two distinct topological sets. One couples 
only to the first and the other to the second order param- 
eter. At the roughening transition only the free energy 
of the cheapest type of steps goes to zero. Simultane- 
ously, the reconstruction order parameter that couples 
to it vanishes, but the other type of order persists inside 
the rough phase. Only the parity type order parameter 
is readily observable in, e.g., x-ray diffraction. There- 
fore, the reconstruction order of the rough phase can go 
unnoticed. For more details, see ref. j^. 

It is conceivable that reconstruction order persists in 
growing surfaces. Imagine a two dimensional (2D) lattice 
with on each site an height variable and a spin degree of 
freedom (representing the reconstruction order). This 
leads to two coupled master equations, one for surface 
growth, e.g., KPZ type dynamics, and another one for the 
reconstruction order, e.g., Glauber type Ising dynamics. 
The local growth probability varies with the Ising config- 
uration and the Ising spin fiip probabilities are affected by 
the local surface height profile. Are these couplings rele- 
vant or irrelevant? If irrelevant, the KPZ sector evolves 
into the stationary non-equilibrium KPZ type state and 
the Ising degrees of freedom reach the Gibbs equilibrium 
state (even while the surface is growing). Coupled mas- 
ter equations of this type have been studied recently in 
the context of specific ID growth models. Those display 
strong coupling between the Ising and roughness degrees 
of freedom, such as growth being pinned down by Ising 
domain walls . We observe a different type of strong 
coupling. 

The 2D restricted solid on solid (RSOS) is one of the 
work horses of surface physics research. Nearest neighbor 
heights differ by at most one, dh ~ 0, ±1. 
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with only nearest neighbor interactions, and dimension- 
less units, K = J/UbT. The J > side of the phase dia- 
gram contains a conventional equilibrium surface rough- 
ening transition ||ic| ] and the growth model version has 
been studied extensively for J > as well P,pl] . 

For J < 0, the model contains one of the simplest ex- 
amples of a reconstructed rough phase , and is prob- 
ably the most compact formulation of the coupling be- 
tween Ising and surface degrees of freedom. The dh ^ ±1 
steps are more favorable than flat dh = segments. At 
zero temperature, the dh = states are frozen out, and 
the model reduces to the so-called body centered solid 
on solid (BCSOS) model. Its surface is rough, but since 
nearest neighbor heights must differ by one, all heights 
on one sublattice must be even and odd at the other, 
or the other way around. This two-fold degeneracy rep- 
resents a checker board type RR phase. The staggered 
magnetization m, defined in terms of the parity spin type 
variables Si — exp^inhj.), is non-zero. 

The dh = excitations appear at T > 0. They form 
closed loops, and behave like Ising type domain walls. 
The reconstruction order m changes sign across such 
loops. Their sizes diverge at the equilibrium deconstruc- 
tion transition Kc = —0.9630 [|lO|. It was found that the 
Ising and roughness variables only couple weakly, i.e., 
that all reconstruction aspects of the transition follow 
conventional Ising critical exponents. Moreover, the ther- 
modynamic singularities in the Ising sector only affect the 
temperature dependence of the surface roughness param- 
eter Kq inside the rough phase. The latter is defined in 
terms of the logarithmic divergence of the height-height 
correlations, 

{{hr+ro-Kor)c^i7rKG)-'Hr). (2) 

The continuum limit confirms this weak coupling. The 
decoupling point of the Gaussian and Ising degrees of 
freedom is there a stable renormalization type fixed 
point 1^. 

We present only our Monte Carlo (MC) simulation re- 
sults in the far from equilibrium limit where evapora- 
tion becomes forbidden. The results look similar closer 
to equilibrium, but crossover scaling phenomena make a 
quantitative analysis more difficult (as expected) . During 
the MC simulation we keep a list of active sites, where 
particles can deposit without violating the RSOS con- 
dition. They are grouped in j = 1, ■ • ■ ,5 sets, accord- 
ing to the only five distinct energy changes AEj that 
can occur during deposition. First we preselect one of 
those 5 sets, with probability {PjNj)/{J2j PjNj), where 

Pj — mm{l,e~^^^) and Nj is the number sites of type 
j. Next, a particle is randomly deposited at one of the 
sites in that specific set j . Rejection free procedures like 
this upset the proper flow of time. We need it, because 
the Metropolis dynamics slows down at low temperatures 
due to an high rejection rate and lack of active sites. To 
restore proper time, we increase the MC time during each 
update step by 1/p x 1/Nj. This reproduces the correct 



value for the KPZ dynamic exponent z ~ 1.6 ± 0.1 at all 
temperatures. 

Fig.|l| shows the susceptibility, x = L'^{{m'^) — (|r7i|)^) 
||l4| , as function of temperature for different system sizes 
Lr. The sharpt maxima seem to confirm the existence 
of a RR phase, but several features are very different 
from equilibrium. The height and width of x/ L'^ do not 
scale with L. At conventional equilibrium transitions, 
the peak height decreases as x/L^ ^t/i^-z ^ 7^-1/4, 
Moreover, the peak position does not converge to a spe- 
cific critical point. Instead it shifts logarithmically, as 
Kp^^y{L) ~ -Aln{L/Lo) with A = 0.77 ± 0.05 and 
Lo = 2.2 ±0.2. 

Next, we monitor the reconstruction order parameter 
TO at _ftr ^ as function of time. It behaves similar as 
in conventional spontaneously ordered phases, but flip- 
flops more frequently than justifiable from finite size ef- 
fects alone. Moreover, the fluctuations in to within each 
phase are too strong. Fig.^ quantifles this in terms of a 
histogram of the number of times a specific value of m 
appears in a typical time series. The distribution has two 
distinct peaks, suggesting the presence of spontaneously 
broken reconstruction order, but the tails have a power 
law shape instead of the mandatory exponential form. 

The above observations suggest the existence of quasi- 
critical reconstructed rough behavior at low tempera- 
tures. The origins of this can be traced to the fol- 
lowing loop dynamics. Fig.(^) can be interpreted as 
a configuration in the ID version of our model, or a 
cross-section of the 2D surface. It shows a domain 
of opposite reconstruction inside an otherwise perfectly 
reconstructed rough configuration. The two fiat seg- 
ments are the domain walls. In equilibrium they move 
apart /towards each other with equal probability because 
deposition/evaporation are equally likely. However, in 
the presence of a growth bias, the defects move more 
likely upwards than downwards. The walls grow in size 
and move up-hill until they get trapped at the top of 
the ridges. Fig.^ shows a trapped loop in an actual 
low temperature MC configuration. Once pinned at the 
ridge line, the loops are slaved to the fluctuations of the 
roughness degrees of freedom. Since the surface fluctu- 
ations are scale invariant (KPZ type in our model), the 
reconstruction order parameter fluctuates critically with 
a power law distribution. Each loop has to follow this 
dance until a new loop nucleates out of the valley and 
annihilates it, or when the encircled terrain happens to 
shrink to zero (flUs-up) by surface growth fluctuations. 

The nucleation of loops takes place in local valleys. 
Consider a deep valley in a perfect BCSOS type RR 
surface conflguration, like in Fig.^ at low temperatures. 
The probability for deposition of one particle at the bot- 
tom of the valley is equal to p = L^^ exp{2K). The next 
event in this local area can destroy the elementary loop 
(by deposition at the same site with probability p = L~^) 
or widen it (by deposition next to it). Annihilation events 
return us to a perfect BCSOS surface that has grown by 
one vertical 2x1 brick. The elementary growth events in 
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the BCSOS model are direct depositions of such bricks. 
In the RSOS model, this process requires an intermediate 
elementary loop excitation state. This implies that the 
time clock in the RSOS model runs slower by a factor 
r = exp(2i4:)(l + 4exp(ii:) + ...) for K < 0. This is the 
origin of the afore mentioned slowing down of the dynam- 
ics for if <C 0. In the following discussion we measure 
time in BCSOS units. 

Loop fluctuations and surface growth events remain 
entangled up to a length scale of about ~ 6 (Fig.(P3)). 
The annihilation of a loop larger than 1^ requires the nu- 
cleation of a distinct new loop from the valley bottom. 
The probability for that is much smaller than for particle 
depositions at the loop itself, which widen the loop and 
make it rise until it becomes trapped on a ridge line. 

The time intervals at which new a macroscopic trapped 
domain of opposite reconstruction order emerges out of 
a valley is independent on the size of the enclosed area. 
Numerically we find t„ ~ exp(— aif), with a — 3.0±0.1, 
independent of loop size. r„ is of the same order of magni- 
tude as simple estimates for the nucleation time of a loop 
of size Ic (r ~ exp(-4if ) in BCSOS time). This part of 
the process is the limiting factor. The second part, in 
which the loop grows into a macroscopic trapped object, 
takes much less time. Positional entropy does not renor- 
malize r„ either. Rough surfaces are scale invariant which 
means that the notion of valley varies with scale. The 
loop enclosed landscape contains many sub-valleys and 
sub-hills, and maybe even an high mountain. However, 
only the deepest valley bottom acts as nucleation site (at 
low temperatures) because loops nucleated in higher sub 
valleys become trapped on sub ridges, and such moun- 
tain lake loops can not grow without additional (rare) 
nucleation events. 

After being trapped on a ridge line, the loop must fol- 
low the growth fluctuations of the surface. Valleys grow 
and shrink (without bias), and fill- up and merge, is the 
life time of a trapped loop of size L on a ridge line subject 
to surface fluctuations only. We expect this time to scale 
as a power law, ~ L^, with z the dynamic exponent 
of the surface roughness degrees of freedom (KPZ like in 
our model). To test this, we measure the decay times of 
large macroscopic defect loops (or order L) as function 
of lattice size L. The data in Fig.(^) collapse indeed on 
one universal curve after rescaling time by Ki . The 
collapse fits best at z = 1.7 ±0.1 (in BCSOS time units), 
which is consistent with the KPZ exponent z = 1.6 ±0.1. 

The above analysis presumes that the nucleation time 
scale T„ ~ exp(— aif) is larger than the surface growth 
time scale, ~ L^. This is valid only well below 
the equilibrium RR transition temperature, and only at 
length scales smaller than, Rc ^ exp(j7f), where loops 
of size L are being nucleated infrequently compared to 
the time scale at which surface growth washes out sur- 
face features of size L. The surface appears as recon- 
structed rough for L < R^. Moreover, the reconstruction 
order parameter appears to be fluctuating in a critical 
manner, since the loops are trapped to the ridge lines. 



and are slaved by the surface fluctuations. So for exam- 
ple, if it were possible to perform x-ray diffraction from a 
growing interface, one would observe not only power law 
shaped peaks associated with the surface roughness. At 
temperatures where Rc becomes larger than the coher- 
ence length of the surface, additional power law shaped 
(critical) reconstruction diffraction peaks will appear. 

At length scales larger than JJc, the surface appears as 
unreconstructed rough. Loops at that large size die by 
nucleation of new loops instead of KPZ surface fluctua- 
tions, and they are not trapped anymore, because loop 
segments of can hop across sub valleys of size R > Rc 
near the ridge line by means of nucleation of new loops 
in those mountain valleys. The peak in the susceptibility, 
see Fig.0, reflects this crossover length Rc- Recall that 
the peak shifts as K — ln(i/Lo). By setting t„ = Tz 
we obtain the same logarithmic behavior. A = z/a is too 
small by about 50%, but this is not a surprise because 
higher order processes renormalize these two time scales 
near Re- 
in conclusion, reconstructed rough phases are absent 
during growth in a strict thermodynamic limit sense, but 
at a more local, and still large length scales (at low tem- 
peratures) the surface grows as if it is reconstructed with 
critical fluctuations in the reconstruction order param- 
eter. Trapping of the loops to the surface degrees of 
freedom at the ridge lines, lies at the core of this. This 
behavior is different from recent results for ID models 
with KPZ and Ising type coupled degrees of freedom. 
There, e.g., the Ising defects become trapped in valleys 
and canyons and pin-down the growth |7|j^. We expect 
to observe crossover to similar structures in 2D by adding 
more interactions in our model and thus vary the local 
growth rates. [|l5|. This research is supported by the Na- 
tional Science Foundation under grant DMR-9985806. 
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FIG. 1. Reconstruction order parameter susceptibility x 
as function of temperature at system sizes L = 8-64. The data 
collapses onto a single curve by the shift K' = K — Kpea.k{L), 
with i4'peak(i) = -0.771n(L/2.2). 





FIG. 3. (a) One dimensional cross section of the surface 
near a valley with two loop segments. The sites a (d) are the 
only active adsorption (desorption) sites. The domain walls 
always move upwards during adsorption, (b) A loop of size 
of Ic nucleated at the bottom of the local valley. Gray and 
white sites have different surface reconstruction parity order. 




FIG. 4. A MC surface configuration (at K = —3.3 and 
L = 64) with a macroscopic loop (the dark line) trapped at a 
surface ridge line. 




FIG. 2. Histogram (shown as insert) of the reconstruction 
order parameter, m, at L = 32 and K = —3.2 from a total 
of 2^* data points using AM = 0.01 as bin width. The tails 



about the peaks at m ■ 
-0.9 ±0.1. 



±1 scale as power laws with exponent 




FIG. 5. Histogram of the decay time of a trapped loop. 
The data collapses onto a single curve by rescaling time by a 
factor L^'^ . 
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